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ABSTRACT 

Aims. In order to improve our understanding of substellar formation, we have performed a compositional and structural study of a 
brown dwarf disk. 

Methods. We present the result of photometric, spectroscopic and imaging observations of 2MASS J04442713+2512164, a young 
brown dwarf (M7.25) member of the Taurus association. Our dataset, combined with results from the literature, provides a complete 
coverage of the spectral energy distribution from the optical to the millimeter including the first photometric measurement of a brown 
dwarf disk at 3.7 mm, and allows us to perform a detailed analysis of the disk properties. 

Results. The target was known to have a disk. High resolution optical spectroscopy shows that it is intensely accreting, and powers 
a jet and an outflow. The disk structure is similar to that observed for more massive TTauri stars. Spectral decomposition models of 
Spitzer/IRS spectra suggest that the mid-infrared emission from the optically thin disk layers is dominated by grains with intermediate 
sizes (1.5 pm). Crystalline silicates are significantly more abondant in the outer part and/or deeper layers of the disk, implying very 
efficient mixing and/or additional annealing processes. Sub-millimeter and millimeter data indicate that most of the disk mass is in 
large grains (>1 mm). 

Key words. Physical data and processes: Radiative transfer, Stars: circumstellar matter, Stars: formation, Stars: individual: 
2MASS J044427 13+25 12164, Stars: brown dwarfs 

I. Introduction 

In the la st decade, severa l studies have confirmed the presence of circumstellar disks around very low-mass stars and brown dwarfs 
(see e.g lApai et alj|200~7l and references therein for a recent review). These objects undergo an accretion phase, like the more 
massive TTauri stars (hereafter TTauri), reinforcing the idea that they mostly form in the same way as stars, i.e. via fragmentation 
and collapse of a molecular cloud. 

Observational evidence for accretion and mass-loss at and below the sub-stellar limit have been obtained for a growing number of 
sources in various star forming regions (e.g [Fernandez & Comeronl 200 1 ; Fe rnandez & Comer6nll2005t iBarrado v Navascues et alj 
2004; iMohantv et ai]|2005t |Whelan et al]|2007l) . iNatta et al.l (120041) have studied the accretion properties of a combined sample of 
20 sub-stellar objects for which the mass-accretion rate is known. These results have shown that the trend of lower accretion rates 
for lower-mass objects continues in the sub- stellar regime. This tr end has been confirmed in subsequent work on larger samples of 
very low mass and sub-stellar objects by e.g Mohant v et al.1 (120051) and lMuzerolle et aL (2001. 

Taking advantage of the unprecedented sensitivity of the Spitzer Space Telescope, it is now possible to investigate th e composi- 
tion of very low mass sta r and brown dwarf disks as previously done for more massive TTauri and Herbig AeBe objects (Apai et all 
20051; iMerin et aLl|2007|) . While the overall disk properties seem to be similar on each side of the sub-stellar limit, significant 
discrepancies are observed when looking at individual targets. 
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In this paper, we present new o bservations o f a you n g brown dwa r f mem b er of the Taurus asso ciation. The mid-IR and millimeter 
excesses reported respectively by [Knap p et al.j (|2004|); iGuieu et all (120071) : IScholz et al.l (|2006) indicate that the object harbors a 
disk. We present and analyze archival sub-millimeter and new mid-IR photometry and spectroscopy data and perform a detailed 
analysis of the structure of the disk via spectral energy distribution (SED) modelling, and of its composition and mineralogy via the 
modelling of the mid-IR spectrum. 

2. Observations 

The target, 2MASS J044427 13+25 12 164, has been initially detected with IRAS bv lKnapp etail d2004l) who identi fied it as a youn g 
stellar object. It was subsequently classified as a M7. 25+0. 25 brown dwarf member of the Taurus association by iLuhmanl d2004). 
In their optical spectrum, they report strong forbidden emission lines indicative of the presence of a jet. The strong Ha emission 
present in their spectrum also classifies the object as an accretor. Now, we have obtained and compiled a comprehensive dataset 
covering the electromagnetic spectrum from 0.43 /urn to 3.7 mm in order to understand the properties of this object. 

2.1. High resolution optical spectroscopy with Keck/HIRES 

We have obtained a high resolution optical spectrum using HIRESr (IVogtet al.lll994h on Keck I on 2006 October 13th. We used a 
setting that covers the wavelength region from 5600 to 10000A, and a slit width of 1.15", yielding a resolution of R=31 000. We 
obtained a single exposure of 40 min. The spectrum was processed and extracted using standard procedures with custom routines 
under IDL (Interactive Data Language). Fig.Q]shows some of the major spectral features observed in the spectrum. 

2.2. Low resolution optical spectroscopy with CAHA/TWIN 

We have obtained low resolution (R=3000, 0.57< A <0.99 fim) optical spectra using TWIN at the Calar Alto Observatory on 2006 
November 20 and 22. We acquired 7 exposures of 10 min each. A photometric standard was observed after the target, but the 
ambient conditions were relatively poor and variable, and the photometric calibration not very accurate. The data were processed 
using standard procedure within IRAF. Fig.|2]shows the average of the seven individual spectra. 

2.3. High angular resolution near-IR images with VLT/NACO 

We used the adaptive optics (AO) NACO instrument dBrandner et al.l2002h on the VLT to image the target at high angular resolution. 
On 2006 December 23 and 24, we obtained images in the H, Ks, NB1.64, NB2.12, and L' filters using the N90C10 dichroic and the 
S13 camera providing a plate-scale of 07013 and a total field of view of »13"xl3". Exposure times were 10 min for each of the 
first four filters, and lh20 min for the L' band. All these observations have been made in jitter mode. The ambient conditions were 
good, with a seeing cr <078 and a coherence time between 10< to <60 ms, acco rding to the ES Q Ambient Condition Databas^B 
The data have been processed using the recommended Eclipse software package (Devillard 1997). A point spread function (PSF) 
reference star (2MASS J04441794+2524513, M4III, Ks=6.0 mag) was observed immediately before or after the target, to be used 
for PSF subtraction. The Strehl ratio in the Ks band reached 16%. 

2.4. VLT/FORS2 imaging 

In order to search for elongations due to the presence of a jet, we have obtained VLT/FORS2 images in the S n filter 
dAppenzeller et al.lll998l) . where the contrast between such extended features and the star is more favorable. The observations 
took place on 2007 February 2 1 . We used the standard resolution 2x2 binning mode leading to a image scale of 0725/pixel resulting 
in a total field of view of 6.8'x6.8'.The pipeline processed image is shown in Fig. [3] The target is unresolved in the T'seeing image. 
No continuum images were taken and no subtraction could be performed. 

2.5. Mid-IR photometry with VLT/VISIR 

We used the mid-infrared imager VISIR dLagage et alj|2004l) at the VLT to obtain images in the [Arm], PAH1, SIV, PAH2, [Nen], 
Ql and Q2 filters. The service mode observations took place between 2006 December 26 and 2007 March 1. Table [4] gives an 
overview of these filters properties. The data were processed with a customary pipeline based on IDL routines. The pipe line includes 
bad p ixel correction and the shift and add of the individual frames, each of them corresponding to a chopping position (Panti n et al.l 
120051) . The target was detected only in the three shortest wavelength (PAH1, Ar in and S iv) filters. The corresponding photometry 
is given in Table|2] The source is not resolved in any of the images. 

2.6. Millimeter and Sub-millimeter observations 

We searched the JCMT public archive to look for for sub-millimeter (sub-mm) data. 2MASS J044427 13+25 12 164 has been ob- 
served with SCUBA on 2004 December 11th at 450 and 850 //m (Program ID: M04BH26A2). We retrieved the data, processed 
them and performed the photometry of the source using standard procedures with the recommended Starlink ORAC software. The 

1 measured in the visible and at the zenith 
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pipeline and its recipes are described in lJenness & EconomouldT999l) . Following the recommendations in the case of weak sources, 
all the bolometers were used and median combined to construct the background. Uncertainties were derived using the "samples" 
method. For each integration, the array moves in a small 3"x3"jiggle grid, in order to account for a a: 1 "offset in the long- and 
short-wave array centers. The "samples" method views each integration as supplying 9 independent data points and retains all 9xN 
data points in its final statistics, where N is the number of integrations (N=160=4 scansx40 integrations in our case). The results 
are given in Table [2] 

2MASS J044427 13+25 12164 was observed on October 12 th , 2007 in time-filler mode at the IRAM/Plateau de Bure interfer- 
ometer (Program ID: rOla, P.I. G. Duvert). The array consisted of 4 antennas in an hybrid compact configuration resulting in a 
beam-size of 9.0x5"oriented North-South, along PA=0°. The receivers were tunes in SSB (LSB) mode at the central frequency of 
86.24GHz (/B.48 mm, corresponding to the J=2-l rotational transition of SiO in its first excited vibrational state), with a 1 GHz 
bandwidth in the two polarizations. Observing conditions were excellent. The flux and error bars reported in Table [2] come from 
a point-source fit in the visibilities at the position of the target. The quality of the observations and the reliability of the IRAM 
pipeline ensure that the point-source fitting measurement does not suffer from any bias. Four quasars with well known fluxes were 
observed for the absolute calibration (0507+179, J0418+380, 0528+134, 3C84). The target is detected, making these observations 
the longest wavelength detection obtained for a sub-stellar object. 

2.7. Mid-IR spectra (5.8-38 (im) with Spitzer/IRS 

2MASSJ04442713+2512164 has been observed with the IRS spectrograph dHouck et al.1 120041) on-board the Spitzer Space 
Telescope in staring mode at both short and long wavelength in the course of program 2 (P.I. Houck) on 2005 M arch 19th. We 
retrieved the public data and processed them using the c2d package which is fully described in Lah uis et~aT1 d2006l) . Two different 
extraction methods are available : full aperture extraction or optimal PSF extraction. The second method is less sensitive to bad 
pixels or bad data samples than the first method providing spectra with less spikes. In this paper we use the spectrum and corre- 
sponding uncertainties obtained with the PSF extraction method. The c2d pipeline furthermore corrects for possible pointing errors 
which can lead to important offsets between different modules. The spectrum ranging from 5 to 38 fim is in good agreement with 
the overlapping VISIR photometry within the uncertainties. The spectrum is shown in Fig. [4] 

2.8. Previously published measurements 

iGuieu et al.l (120071) presented IR AC and MIPS photo metry of 2MASSJ04442713+2512164. It was the only brown dwarf of their 
sample with a detection at 70 ^m. lScholz et al.l (12006) presented 1.3 mm observations of 2MASSJ044427 13+25 12 164, and derived 
a disk mass significantly larger than any other disk mass in their sample. The corresponding photometric measurements are quoted 
in Table [2] 

3. Stellar parameters 

Table|5]gives a summary of the stellar parameters of 2MASS J044427 13+25 12 164. 

Spectral type. - Using the low resolution optical spectra obtained with TWIN, and the classification scheme described in 
Mart in et alj ([1999), we derive a PC3 index of 1 .64, corresponding to a spectral class of M7. 1 . Fig.[2]s hows indeed that the spectrum 
is very similar to that of the field ultra-cool dwarf vB8, which has been classified as an dM7 (Marti net alJI 1999h . The differences 
between vB8 and 2MASS J044427 13+25 12 164 spectra can be ascribed mainly to gr avity effects. We thus attribute a final spectral 
type of M7. 0±0.5. These results are in good agreement with the previous estimates of lLuhmanl (|2004), who reported a spectral type 
of M7.25+0.25 and a T eff =2838 K. 

Rotation velocity. - We derive the rotation velocity (v sin i) of the target by cross-correlating the HIRES spectrum with a g rid of 
"spun-up" templates of a slowly rotating standard of similar spectral type obtained the same night (vB 10 Basri & Reinersl2 006). The 
cross-correlation was restricted to orders free of strong stellar emission lines, strong telluric features, or gravity-sensitive features. 
The best match is obtained for v sin i =12+2 km s . 

Veiling and exti nction. - From our optical spectrum we derive an extinction Ay=0, consistent with the previous estimates Aj=0 of 
iLuhmanl d2004l) . We also estimate the optical continuum veiling r = r F 6750 \/'F' 6750 Xcont> ^ ue *° excess emission from the accretion 
shock, by comparing the strength of the A615Q A TiO molecular band in our high-resolution spectrum to that of a template of 
similar spectral type and luminosity type (vBlO). The best match is obtained for a veiling of 0.4+0.1. 

Accretion and mass loss. - Both the low resolution CAHA/TWIN and high resolution Keck/HIRES spectra show strong emis- 
sion lines superimposed on the stellar continuum (Fig. Q] and [2] and Table [3]). These lines are characteristic of intense accretion 
process: the H a and Can triplet lines show equ i valent widths typical of very active classical TTauri stars. A strong level of H a 
emission had already been reported by ILuhmanl d2004l) but no value of equivalent width was reported in their paper. The H„ line 
has a broad and asymmetric profile, as shown in Fig.Q] We tentatively derive a mass accretion rate using the method described in 
Mohantv et alj d2005|) base d on the pseudo-equivalent width of the Ca n emission at /18662 A, and using the spectral type and T e ff 
reported in ILuhmanl (120041) (M7.25, T e fj =2838 K) and the corresponding mass for an age of 5 Myr as given by the DUSTY models 
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Table 1 Explored parameters in our grid of models and results of the Bayesian analysis. 



Parameter 


Min. value 


Max. value 


1 ' sample 


sampling 


Valid range 


R r AT Tl 


JUU 




1 

1 


fixed 




Mdust LMqJ 


1U 


rn-5 
1U 


i 
1 


fixed 




a max [mm] 


1 


1 


1 


fixed 




inclination [°] 





90 


10 


linear in cos ( 


< 55° 


surface density exp. 


-2 





5 


linear 


-2 to -0.5 


flaring index 


1.0 


1.25 


6 


linear 


1.1-1.2 


R in [AU] 


0.015 


0.15 


6 


log 


0.02-0.1 


Ho [AU] 


20 


80 


13 


linear 


30-60 


stratification 


0. 


0.1 


5 


linear 


no constraint 



(M=0.045 Mgl lBaraffe et~ai1l2002l) . W e obtain a value of logM=-9.7 Mp^ yr of t he same order than observed for very low mass 
stellar and sub-stellar objects in Taurus dMuzerolle et aHl2003tlWhite & B asri 2003) and other associations of comparable ages (e.g 
Barrado v Navascues et al. 2004; Moh antv et alj|2005l) . We note that the Can emission lines are observed in both HH objects and 
strongly accreting stars, so that the fl ux of this line might be common to both accretion and mass-loss and the mass accretion rate 
derived corresponds to an upper limit dAlencar et al.ll2.Q05b . At the same time, numerous, strong forbidden optical lines (O i, N n and 
S n) suggest that the object is also undergoing mass loss. 



4. Disk SED modelling 

We compare the SED of the ta rget with a grid of 45 000 passive disk models computed with a Monte Carlo Radiative Transfer code 
(MCFOST), fully described in Pi nte et aTT ( l2006l) . In short, the code computes the disk temperature structure by propagating photon 
packets originally emitted by the central star through a combination of scattering following Mie theory, absorption and re-emission 
until they exit the computation grid and reach the observer. After the temperature computation, SEDs of the disk are estimated by 
emitting and propagating the proper amount of stellar and disk thermal photon packets, with a weight ensuring that each wavelength 
bin is properly converged. Ensuring energy conservation is a key issue for radiative transfer codes because of the very large surface 
densities reached in sub-stellar disk inner regions. Our two-steps calculations output two SEDs calculated by generating packets in 
two in dependent ways. The accuracy has been systematically a posteriori validated by comparing these two SEDs (see lPinte et al.l 
2006, section 3.3 for more details). Synthetic SEDs are computed simultaneously at all inclinations in 10 bins regularly spaced in 
cos i, with 60 wavelength bins logarithmically spaced between 0. 1 and 4 000 yum. The temperature computation is performed using 
6.4xl0 5 packets, whereas synthetic SEDs are computed using approximately 6.4xl0 6 photons. 

The parameter space for the disk model is strongly degenerated making it impracticable and too much time-consuming to 
explore a grid of sufficient resolution over all parameters. In order to reduce the size of the parameter space to be explored, we run 
preliminary models to get c onstraints on some of the parameters. The general shape of the silicates features is roughly reproduced 
by amorphous olivine dust dDorschner et al.|[l995l) . In the following global SED modeling approach, we do not try to get a perfect 
match of these features and adopt a composition of 100% amorphous olivine. A detailed analysis of the mineralogy is performed 
afterwards. 

As input we use the stellar parameters (T e ff=2838 K, Ay=0.0 mag, L=0.028 Lq) d erived bvlLuhmanl d2004l) . the corresponding 
radius (R=0.64 R©), the photospheric emission predicted by the NextGen modelfl of Bar affe et al.l (120021) . and an assumed outer 
disk radius of 300 AU. We assume a grain size distribution following a power-law with dn{a) oc a~ 3 ' 5 da, between a m ; n =0.03 //m 
(chosen very small so that its exact value has no effect on the observables) and the maximum grain size a max considered as a free 
parameter. The spectral index between 1.3 and 3.47 mm log(F v )/ log(v) = 2.7 indicates the presence of grains significantly larger 
than interstellar ones. A maximum grain size around 1 mm gives a good fit of the spectral index and we adopt a max = 1 mm. With 
this dust population, the millimeter emission is well reproduced with a dust mass of 10~ 5 M Q . Assuming a typical dust to mass 
ratio of 100, this le ads to a total disk mass of 10~ 3 M in good agreement with the previously published estimate of 2.15xl0~ 3 Mq 
dScholz et alJl2006l) . 

The star properties, disk outer radius, grain size distribution, dust composition and dust disk mass are kept fixed in the rest of 
the modelling. We then run a grid of models that systematically samples the remaining parameters: disk inclination, inner radius, 
scale height, flaring and surface density exponents. The explored ranges for each parameter are described in TableQ] 

We performed a comparison of the synthetic SEDs with observations using a y 2 minimization. The best model is presented in 
Fig. |U It gives an overall good match of all observational points from optical to millimeter regimes. As already mentioned, we 
did not intent to get a quantitative fit of the silicates emission bands but rather focus on the strength of these features. A Bayesian 
analysis is performed using the reduced x 1 to get an estimate of the validity range for each of the parameter explored dPress et al.l 
1992; La v et al.ll 1997h IPinte et alj|2007l) . A relative probability exp(-^f 2 /2) is calculated for each individual model. The relative 
likelihood for each of the parameters is obtained by adding the individual probabilities of all the models with a given parameter 
value. These probabilities are the result of a marginalization of the parameter space successively over all dimensions and not a cut 
through the parameter space. In that sense, they account for potential correlations and interplay between parameters, and quantitative 
error bars can be extracted from them. The results are presented in Fig. [6] 



2 for T cff =2 800 K 



H. Bouy et al.: Structural and compositional properties of brown dwarf disks: the case of 2MASS J04442713+2512164 



5 



5. Disk mineralogy 

The features present in the Spitzer IRS spectrum probe emission originating essentially from silicate dust grains less than a few 
microns in size (larger grains are featureless) and located in the optically thin surface layers of the inner disk (<2 AU). The following 
section presents a detailed analysis of the properties of these grains. 

The data show several evidences of grain growth and processing. The slope of the SED at longer wavelengths indicates that large 
grains («1 mm) must be present in the cold disk, while the trapezoidal silicate emission profiles seen in the Spitzer IRS spectrum 
around 10 and 18 fim are typi cal of grain growth to micron size and/or pres ence of crystalline silicate forsterite, the Mg-rich end 
member of the olivine group (iBouwman et alj|200lt van Boe kel et al.l l2005). The IRS spectrum shows clear emission features at 
around 11.3, 23-24, 27-28, 33 microns that we attribute to crystalline silicate grains, in particular crystalline forsterite which shows 
features at 10.1, 11.3, 16.2, 19.7, 23.8, 27.6 and 33.6 fim, and crystalline enstatite (Mg-rich pyroxene) which shows features at 
9.3, 10.6, 1 1.7, 19.6, 21.6, 23.0, 24.5 and 28.2 fim. The presence of crystalline grai ns in emission indic ates advanced processing of 
the circumstellar material in the disk compared to the diffuse interstellar medium (Kemper et al. 2004). In order to investigate the 
composition, distribution and level of processing of the grains of the inner disk (inner meaning here an d after ^2 AU) of t he target, 
we have followed the procedure used in interpreting the Spitzer IRS spectrum of a very low mass star bv lMerfn et al.l (120071) . Details 
of the procedure and references can be found therein. Briefly, the method includes three major steps: 

1 . subtraction of a continuum 

2. fit of the continuum-subtracted 10 fim amorphous silicate feature (between 7-15 fim) using a grid of pre-calculated absorption 
efficiencies for five types of silicate species assumed to emit at the same T e fr ■ A x 1 minimization provides the relative silicate 
mass fractions and best blackbody T?jj?. 

3. fit of the residual spectrum (between 15 and 35/mi) as in 2., providing a second, cooler blackbody T™ w and relative mass 
fractions for the cold component. 

Following Merin e t al.l d2007l) and references therein, the models includes five major grain species identified in circumstellar 
disks (3 types of amorphous silicates: olivine, pyroxene and silica; and 2 kinds of crystalline silicates: forsterite and enstatite), and 
assumes two representative (from a spectroscopic point of view) grain sizes (0.1 and 1.5 fim). The contribution of larger grains 
cannot be studied usi ng the curren t dataset and would have been fitted in the continuum. 

As mentioned in iMerfn etaTI d2007l) . the main limitation of this analysis lies in the degeneracy of the contributions of the 
grain species. The silicate emission indeed contributes to the emission over almost the entire IRS spectrum, preventing to measure 
accurately the purely photospheric+disk continuum emission which we need to subtract in step 1. To assess the uncertainties and 
limitations of this analysis, we investigate the mineralogy of the disk using independent methods for the continuum normalization. 

5.1. Continuum normalization 

MCFOST disk model continuum. - Our first approach consists in using the best-fit SED obtained with the MCFOST code as the 
continuum. The main caveat is that a treatment of the optical properties of amorphous olivine is included in the MCFOST models 
(see Fig. [4] the 10 and 20 fim amorphous silicate features are relatively well fitted by the model), preventing to use the disk model 
as it is. In order to make a more reliable estimate of the continuum, we have computed the MCFOST models without any treatment 
of the amorphous olivine emission by replacing the optical constants of this species with a power-law interpolation between 5 and 
30 fim. The 10 and 20 fim silicate features are removed, but the energy originally re-emitted in the silicate features is now spread over 
the spectrum (the disk absorption indeed remains the same) and affects the continuum, as shown in Fig. [4] In order to compensate 
this effect, we estimate the continuum using a smoothed version of this synthetic SED, shifted to match the 10 fim feature. The 
corresponding mineralogical analysis gives an overall good fit of the different features (see Tables [6] and 13 . The best fit is obtained 
for a mass fraction in crystalline grains of 25%. Fig.|7]gives an overview of the results. 

Polynomial fit. - We also d erive an estimate of the continuum using a Lagrange polynomial fit following the method described in 
iKessler-Silacci et all d2006l) (in their Section 4.1, case 1). The fit obtained using this continuum gives a mass fraction in crystalline 
grains of 31%, consistent with the previous estimate. The corresponding mineralogical analysis gives an overall good fit of the 
different features as shown in Fig. [7] (right panels) and by the reduced-^ 2 reported in Tables [6] and [7] As an additional test, we 
computed a second polynomial continuum with slightly different properties, leaving more flux around 35 fim. All the features, even 
the one at 33 fim, are well fitted assuming this second continuum. The estimated cristallinity, «20%, is significantly lower than that 
estimated with the first polynomial continuum. Leaving more flux at longer wavelength naturally favors amorphous silicates, leading 
to a lower crystallinity fraction. The spectral features associated with amorphous silicates are indeed broader than that associated 
with crystalline silicates, and are therefore better suited to fill the flux deficit at long wavelengths in this continuum. 

While the first polynomial fit of the continuum gives overall results very similar to those obtained with the MCFOST models 
(see Fig.[7j>, the second polynomial fit gives a similarly good fit but with very different results, illustrating the difficulty of choosing 
the right continuum. Fortunately, the synthetic SED obtained with MCFOST provides a first guess of the shape of the continuum, 
or at least of plausible range of acceptable continuum. The second polynomial fit is significantly different from the synthetic SED 
obtained with MCFOST, and therefore less likely to be representative of the true continuum. It nevertheless provides an interesting 
check and a lower limit on the crystallinity level, at about 20%. This iterative analysis, which takes advantage of both the MCFOST 
and polynomial fits, leads us to conclude that the crystallinity mass fraction cannot reasonably be larger than ^35% or lower than 
^20%, with the best fits obtained in the range between 25-31%. The results of the analysis are summarized in Tables [6] and [7] 
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5.2. Grain sizes in the surface layers of the inner disk (<2 AU) 

Interestingly, the best fit is obtained for a mixture made of essentially intermediate size grains. We investigate the effect of the grain 
sizes on the results of the fit. The method originally uses two grain sizes: small grains a sma n=0.1 fim, and intermediate size grains 
<2big=l-5 i-im. Increasing the maximum grain size to 3.0 fim gives a better fit of both the blue side of the 10 and 20 fim features, but a 
worse fit of their red sides, of the 29.2 and 30.9 fim features and of the relatively strong crystalline forsterite 33 fim feature. Larger 
grains (>6 fim) did not lead to any acceptable fit. We tentatively used 3 grain sizes, but the number of free parameters becomes too 
large and the results inconclusive. I n conclusion, both the high a nd low temperature components are best fitted with au g = 1.5 fim 
grains, consistent with the results of IKessler-Silacci et al. (2006) for more massive TTauri M stars. Smaller a^g produces narrower 
features, and larger a^g a poor fit of the red side of the different features. For the rest of the mineralogical analysis, we adopt a 
maximum silicate grain size of at,ig=l-5 fim. Finally, we stress that the 27.0 fim feature was never well fitted by any of the models 
mentioned above, must be associated to species not included in our analysis. 

6. Discussion 

6.1. Disk structural properties 

The unique coverage and sampling of the SED allows to constrain a number of disk structural properties. 
2MASS J04442713+2512164 disk is th e most massive dis k known to date around a sub-stellar object, surpassing other observed 
brown dwarf disks by a factor 4 or more (IScholz et alJ20 06). This recalls the results observed for more massive TTauri stars showing 
that circu mstellar disks of o bjects with well developed outflows are significantly more massive than disks of sources with no such 
outflows dCabrit et al.lll990l) . Although no outflows was resolved in our NACO images, its presence is suggested by the numerous 
strong forbidden lines present in the optical spectrum, as mentioned above. 

Fig. [6]presents the relative figure of merit estimated from the Bayesian inference of each disk parameter probed by our SED 
analysis with the MCFOST models. We obtain quantitative constraints on the inner radius, scale height and flaring exponent of 
the disk. The inner radius is constrained between 0.02 and 0.1 AU with a peak probability around 0.04 AU. The flaring index is 
constrained between 1 . 1 and 1 .2 and flat models j3 = 1 are excluded. Acceptable values for the scale height are obtained in the range 
30-60 AU, with a most probable v alue close to 45 AU . This large value for the scale height is consistent with the hydrostatic scale 
height expected for a brown dwarf dWalker et alj|2.004b . The best model gives indeed a mid-plane temperature of 12.5 K at 100 AU. 
Assuming a central star mass of 0.045 Mq and a vertically isothermal disk (which is not the case in the model), the mid-plane 

temperature translates in a scale height ^kB^Ti^/GM^fi as 35 AU. The increasing temperature towards the disk surface will result 
in slightly a larger "effective" scale height, therefore in good agreement with the scale height deduced from the SED fitting. All 
these properties are very similar to the properties of disks of more massive TTauri, suggesting that similar processes are at work on 
each side of the sub-stellar limit. 

Some of the explored parameters cannot be constrained by our modelling of the SED, resulting in flat probabilities. The settling 
level is almost unconstrained and good fits can be found with and without settling. This is consistent with the fact that the dust 
scale height is in agreement with the gas scale height estimated assuming hydrostatic equilibrium and indicates that the dust and 
gas are probably well mixed within the disk. The surface density exponent is also only loosely constrained. Models with flat surface 
densities are excluded, the opacity in the central regions of the disk becoming too low to produce enough near-infrared excess. The 
fit is also insensitive to the disk inclination as long as the star remains directly visible by the observer. When the inclination becomes 
too large (> 55°), the star starts to be occulted by the disk and the optical and near-infrared fluxes are strongly reduced. 

We also tentatively investigate the external disk truncation via the influence of the outer radius on the fit by computing synthetic 
SEDs for three different values of R out . Fig. [5] shows synthetic SEDs obtained for R out =300, 10 and 1 AU. The SED remains 
almost unchanged down to an outer radius as small as 10 AU, while keeping the disk mass constant. At R out =1 AU, the flux in 
the sub-mm/mm range decreases dramatically, the disk becoming optically thick, and indicating that in this case the limit must lie 
between 1 and 10 AU. This emphasizes a limitation of such analysis and the fact that the outer radius and disk truncation cannot 
be directly constraine d from the sub-millim eter/millimeter photometry. A way to break this degeneracy is to image the disk, as 
illustrated recently in Luhm anet al.l (120071) for another sub-stellar member of the Taurus association. We obtained deep adaptive 
optics images with the aim to search for extended diffuse emission coming from the disk or the jet. The faintness of the source, and 
the intermediate inclination of its disk made these observations and the data analysis challenging. The current results are negative 
and no obvious extension is found in the images. 

6.2. Crystallization of Silicates 

Table [6] gives a summary of the main properties of the 0.1< a <1.5 fim population of grains in the hot and cold temperature 
components of the optically thin layers of the inner disk (^2 AU) probed by the mid-IR Spitzer spectrum. Both components seem 
domi nated by intermediate size grains {a— 1 .5 fim), indicating a relatively high level of processing. This is consistent with the results 
of IKessler-Sil acci et al.l (120061) who report a statistical dependency of grain size on spectral type for more massive TTauri stars, 
cooler objects harboring disks made of larger grains. 

The hot and cold temperature components probed by the IRS spectrum are produced in different locations of the inner disk, the 
cold component corresponding mostly to the outer part and/or deeper layers of the inner disk, and the hot component corresponding 
to the surface layers and/or inner region of the inner disk (see Fig. [8]). Interestingly, the cold temperature component contains 
a significantly higher fraction of crystalline species than the hot component (10 to 15 times larger). Since the cold component 
contributes to about half of the mass, as shown in Table [7] it implies that crystalline species are mostly located in the deeper 
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layers/outer region of the disk. This result does not fit in the classical picture where thermal annealing of amorphous silicates takes 
place by means of heati ng to high temperatures very close t o the central star (at about 8 00 K, corresponding to 10 _1 ~10~ 3 AU in the 
case of brown dwarfs, iGailll 1998c lHallenbeck et al.ll2000t iKessler-Silacci et alj|2007l) and resulting in a higher crystallinity closer 
to the star than in the outer parts of the disks (Ivan Boekel et aU2004l) . In the following sections we discuss our results in the context 
of tree scenarios that may be responsible for the production of crystalline silicate grains in protoplanetary disks: (1) annealing in the 
very early stages of the accretion disk formation; (2) annealing in the hot inner regions of the disk followed by large-scale outward 
radial transport; (3) local annealing in the outer cold regions produced by energetic events such as flares or shocks. 

6.2.1 . Thermal annealing in the early phase of disk formation 

iDullemond et alj ([2006) suggest that dust crystallization events happen in the very early stages of the disk formation and relate the 
abundance of crystalline species in the entire disk to the rotational velocity of the parent cloud. Rapidly rotating clouds produce 
massive disks with lower accretion rates and lower crystallinity, while slow rotating clouds produce less massive disks with lower 
accretion rate but high levels of crystallinity. This process would fill the protoplanetary disk with a relatively homogeneous radial 
distribution of crystalline species. Fig. [10] shows several cases that do not fit in this scheme (larger accretion rates but significant 
levels of crystallinity), including 2MASS J04442713+2512164 which additionally has the most massive disk known to date around 
a brown dwarf. If at work, this scenario must therefore be followed by additional annealing and/or mixing processes as the disk 
evolves with time, and the question arises as to what explains the larger fraction of crystalline species among the 0.1< a <1.5 jum 
population of grains in the deeper layers/outer regions of the inner disk probed by our study. 

6.2.2. Thermal heating in the disk inner regions and radial mixing 

Boc kelee-Morvan et al.l d2002l) and then lKeller & Gaill d2004l) described models were turbulent diffusion carries crystalline silicates 
from the inner to the outer disk, while the surface layers are slowly spiraling inwards to reach high temperatures and anneal. 
The accretion rates required in their simulations to transport the annealed grains to the outer regions were inconsistent with the 
most recent models of formation of the solar system. Slower accretion rates were not able to reproduce the observed abundances 
of processed grains at large radii mostl y becau s e the outward transport of processed grains was surpassed by the inward flow of 
material accreting onto the central star. ICieslal d2007l) recently developed further these studies and presented a two-dimensional 
protoplanetary disk model including the treatment of both the radial and the vertical transport of grains. The new simulations show 
that the outward transport of processed grains is a stratified mechanism which appears to be more efficient than initially predicted by 
the previous one-dimensional models. Outward transport occurs most efficiently near the disk mid-plane, where the pressure due to 
the gradient of density eases the transport to the outer regions, and where the grains do not have to fight the inward flow of accreting 
material. Accretion plays a key role in these scenarios. The efficiency and spatial scales on which these additional mechanisms are 
taking place are difficult to quantify with the current data, accretion being a strongly variable phenomenon, and the scales probed 
by mid-IR spectroscopy depending on the geomet ry of the disk (see discussio n in section |6"3l l. However, recent statistical studies 
carried over a large sample of solar mass TTauri by Kessl er-Silacci et all (120071) showed no correlation between the grain properties 
and the accretion as measured by the H a equivalent width. Fig. [K)]shows the crystalline mass fraction of samples of Herbig Ae/Be 
stars, TTauri stars and brown dwarfs found in the literature. Even though the limited size and the inhomogeneity of the samples 
(spanning different ages, star forming regions and measurement methods) represented in this figure prevent us to draw any firm 
conclusion, there is no indication of a correlation between the crystalline mass fraction and the accretion rate on either side of the 
sub-stellar limit. 

Finally, if the mixing is not efficient enough, the lower crystalline mass fraction in the hot component of 
2MASS J044427 13+25 12 164 could be explained by reverse transformation (melting) of crystalline to amorphous species if the 
newly formed crystals are not transported fast enough to deeper layers or outer radii. 

6.2.3. Transient grain annealing 

Energetic events (shocks in the disk or flaring of the central object) could be responsible for a s ignificant produc t ion of crystalline 
silicates in outer regions where the temperature is much cooler than the annealing temperature. Harker & Deschl d2002l) computed 
the elevation of temperature due to shock waves triggered by gravitational instabilities in the case of typical TTauri systems and 
concluded that annealing of silicate grains in the 5-10 AU region is possible. Because of their smaller thermal resistance, small 
grains are more likely to be annealed and the fraction of crystalline species is therefore expected to be larger for small grains that for 
large grains. Our analysis of the mid-IR spectrum of 2MASS J044427 13+25 12 164 shows that the small grains of the hot components 
are indeed mostly crystalline, but that they are mostly am orphous in the cold component. Spectroscopy at longer wavelengths (with 
e.g the Herschel Space Observatory, Pilbr att et al.l 12004) should bring new insights on the contribution of transient annealing in 
protoplanetary disks. 

6.3. Additional crystallization processes 

Mols ter et al.l d 19991) have discussed several mechanisms to increase the degree of crystallinity of red giant disks without the need 
for radial mixing: selective removal of amorphous grains (either by transport or coagulation) or additional crystallization methods 
(below the annealing temperature). Their conclusion is that the latter explanation is the most likely one, selective transport occurring 
on much slower timescales and coagulation affecting crystalline and amorphous species in the same manner. Because of the simi- 
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larities between red giant and young star disks, they also conclude that such additional mechanism(s) could be at work in young star 
disks. Recent laboratory experiments have shown that exothermic chemical reactions involving the graphitizati on of a carboneou s 
mantle layer covering amorphous silicate grains are able to transform these grains into crystalline silicate grains dKaito et al]|2007l) . 
Both such organic species a nd crystalline silicate grains have been found in the comet 8 IP/Wild 2 samples retur ned by the Stardust 
mission feller et al.l l2006) and in the observations of comet 9P/Tempel 1 after the Deep Impact experiment dHarker et al.| [2007; 
iTozzi et alj|2007l) . 



6.4. Grain growth 

Coagulation into large millimeter size grains is clearly demonstrated by the sub-millimeter and millimeter fluxes of 
2MASS J044427 13+25 12 164. Millimeter photometry mostly traces the thermal emission of cold grains in the outer part (>20 AU) 
of the disk mid-plane, as shown in Fig. [8] The presence of silicate features in the mid-IR spectrum shows that intermediate size 
grains (a «1 /urn) are also still present in the inner disk (<2 AU, see Fig . [HJ . On the other hand, the same mid-IR spectral analysis 
shows a lack of small grains (a «0.1 /urn) in these regions/surface layers (see Table [6]). These three results together with the much 
higher level of crystallinity found in the cold component of the inner disk show that the disk of 2MASS J044427 13+25 12 164 has 
reached a relatively advanced level of processing, leading to the formation of elaborate crystalline structures and large grains. 

Unfortunately, the non-detection of the disk in direct imaging and the degeneracy of the SED modelling of the optically thin 
disk at sub-mm and mm wavelengths do not allow to derive any constraints on the vertical and radial distribution of the large grains. 
Fig. [6] shows that in the current state of the data and of the models, no constraints can be derived on the vertical distribution of 
the grains. Simple considerations on the dynamics of large dust grains, which is mainly determined by gravitational sedimentation 
towards the disk mid-plane, naturally sugg est a stratified disk structure with large grains mostly located in the disk mid-plane, and 
smaller grains at the surface of the disk (Chiang et al. 2001). This simple picture can nevertheless be complicated by radial and 
vertical mixing, and fragmentation of the grains via collisions. The larger crystalline mass fraction of the cold component suggest 
that such settling is at work and must be very efficient. Direct imaging of the disk at different wavelengths and/or interferometric 
observations are required to constrain the grain size radial and vertical distributions. 



6.5. Comparison with previous studies and limitations of the analysis 

The methods used to model the dust properties suffer from several limitations. None of the models in this paper (or in the literature) 
provide a self-consistent and physically correct dust model, mainly due to several unconstrained parameters or to simplifications. 
These include the limited number of temperatures (two), the limited number of grain species (both in terms of sizes and nature) 
and in particular the important effects of dust species without infrared features on the continuum, cold grains, and the optically thin 
continuum from large (>10 fim) silicate grains that do not show up as spectral features, and the porosity of the grains. The values 
derived in this paper should therefore be taken with great caution. Keeping these limitations in mind, our analysis nevertheless 
allows to perform qualitative comparisons with other objects observed and studied using similar methodology. 

Fig. [9] compare s the mass fraction i n crystalline grains of 2MASS J044427 13+25 12 164 to that of mo re massive 
Herbig AeBe s tars (Ivan Boekel et all 1 20051) . TTauri stars dMeeus et ail 120031; ISargent et all l2006t iMerin et alj 120071) and other 
brown dwarfs (Ap ai et al.l 120051) . With a crystalline mass fraction between 20-30% and a S 11.3/S9.3 ratio equal to 1.05, 
2MASS J044427 13+25 12 164 has an overall degree of crystallization comparable to objects of similar mass and ages. 

The compositional analysis of the M5.5 star Par-Lup-3-1 performed by Meri n et ail d2007l) provides an opportunity to compare 
the results obtained for two objects close to the sub-stellar bounda ry and with similar ages (the estimated ages of the Lupus III 
and Taurus associations being around 1-3 Mvr). lMerfn et al.l ([2007) report a much larger fraction of small grains (0.1 fim) in both 
the hot and cold components, contrasting with the lack of such small grains found in our analysis of 2MASS J044427 13+25 12 164 
spectrum. The millimeter observations being not available for Par-Lup-3-1, the following comparison between the properties of the 
two disks holds for the inner disk (<2 AU) probed by the mid-IR analysis. The hot component of Par-Lup-3-1 contains most of the 
crystalline mass, while we find a much higher crystalline mass fraction in the cold component of 2MASS J044427 13+25 12 164. The 
most striking difference b etween these two objects lies in the lack of any accretion signature in the optical spectrum of Par-Lup- 
3-1 dComeron et al.ll2003l) . while 2MASS J044427 13+25 12164 displays one of the strongest accretion related Ha luminosities. 
Fig. [10] shows nevertheless that there seem to be no correlation between the accretion luminosity and the crystallinity of sub- 
stellar and very low mass objects. Several other independent parameters could well be responsible for the differences observed 
in 2MASS J044427 13+25 12 164 and Par-Lup-3-1 disk compositions. The initial environmental c onditions or the disk p rocessing 
mechanisms at work on each side of the sub-stellar limit must indeed play crucial roles as well (Dulle mond et al.l [20061) and are 
not taken into account in the current analysis. The disk geometries, and in particular the orientation, must also affect the results of 
our analysis. The smaller inclination of 2MASS J044427 13+25 12164 (i <55° instead of i >70° for Par-Lup-3-1) possibly implies 
that the mid-IR spectrum probes larger depths into the low-density surface layers, "artificially increasing the overall mass fraction 
in the cold component when compared to Par-Lup-3-1. On the other hand, a smaller inclination also corresponds to a smaller 
projected surface of the internal edge, where most of the mid-IR emission probed by the mid-IR Spitzer spectrum originates. 
These opposing geometrical effects are not fully taken into account in the IRS spectral analysis, and not only complicate a direct 
comparison between different sources but also add an unknown weighting to the overall compositional analysis. Similar analysis 
on large samples covering a statistically well defined range of geometries are required to draw more quantitative conclusions and 
comparisons. 
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7. Conclusions 

We have presented a complete analysis of a brown dwarf SED from the optical to the millimeter, and a detailed compositional 
analysis of its entire (5-35 fim) mid-IR Spitzer spectrum. The conclusions are severalfold: 

1. The disk has geometrical properties similar to that observed for more massive TTauri. In the current state of the models and of 
the observations, the analysis of the SED does not allow to derive any constraints regarding the outer radius of brown dwarf 
disks. Millimeter observations lack the angular resolution to study disk truncation. Additional observations including direct 
imaging, and spectroscopy of the inner region via interferometry or e.g with the future Herschel Observatory, are required in 
order to study the vertical and radial distribution of grains, but also to search for evidence of disk truncation. 

2. The sub-mm and mm photometry indicates that most of the mass of the disk is in grains significantly larger (<1 mm) than in the 
diffuse interstellar medium (100A-0.2 jum). Grain growth to millimeter size grains can therefore occur as rapidly as 1-3 Myrs. 

3. As for the very low mass M5.5 star Par-Lup-3-1, the mid-IR Spitzer spectrum is better fitted using two temperature components 
rather than a single one. The hot component displays a significantly lower crystallinity than the cold component, indicating 
an advanced level of processing, and possibly contrasting with Par-Lup-3-1. This also suggests that: (i) the radial and vertical 
mixing of crystalline species in the disk could be very efficient; (ii) that crystallization processes other than high temperature 
annealing close to the central object might be at work; (iii) that crystals might be melted back to amorphous state faster than 
they are transported out of the hot layers where they form 

4. In the current state of the observations, the composition of very low mass stellar and sub-stellar disks does not seem to correlate 
with accretion. 

5. The proportions of intermediate size grains (1.5 //m) in both the hot and cold components are larger than those derive for Par- 
Lup-3-1, providing additional evidence of a more advanced level of processing than the M5.5 very low mass star even though 
the two objects are located in associations of similar ages (1-3 Myrs). 

These results provide additional evidences that brown dwarfs disks are similar to those of more massive TTauri stars, and that 
the initial steps of planet formation are present during the process ( Apai et al. 20051). 

Acknowledgements. H. Bouy acknowledges the funding from the European Commission's Sixth Framework Program as a Marie Curie Outgoing International 
Fellow (MOIF-CT-2005-8389). We are grateful to our referee Daniel Apai for his comments and advice which helped improve this manuscript considerably. 

This work is based on observations obtained at the VLT which is operated by the European Southern Observatory. This work makes use of DENIS data. The 
DENIS project has been partly funded by the SCIENCE and the HCM plans of the European Commission under grants CT920791 and CT940627. It is supported 
by INSU, MEN and CNRS in France, by the State of Baden- Wiirttemberg in Germany, by DGICYT in Spain, by CNR in Italy, by FFwFBWF in Austria, by 
FAPESP in Brazil, by OTKA grants F-4239 and F-013990 in Hungary, and by the ESO C&EE grant A-04-046. This publication makes use of data products from 
the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute 
of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This work is based [in part] on archival data 
obtained with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. The 
IRS was a collaborative venture between Cornell University and Ball Aerospace Corporation funded by NASA through the Jet Propulsion Laboratory and Ames 
Research Center. The authors acknowledge the data analysis facilities provided by the Starlink Project which was run by C CLRC on behalf of PPARC. This work has 
made use of the Vizier Service provided by the Centre de Donnees Astronomiques de Strasbourg, France lOchsenbein et al. 2000). This research used the facilities 
of the Canadian Astronomy Data Centre operated by the National Research Council of Canada with the support of the Canadian Space Agency. Some of The 
data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the 
University of California and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W.M. 
Keck Foundation. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Mauna Kea has always had 
within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. Some of The data presented 
herein were obtained at the Calar Alto Observatory, which is operated jointly by the Max-Planck-Institut fur Astronomie (MPIA) in Heidelberg, Germany, and the 
Instituto de Astrofi'sica de Andaluci'a (CSIC) in Granada/Spain. 



References 

Alencar, S. H. P., Basri, G, Hartmann, L., & Calvet, N. 2005, A&A, 440, 595 
Apai, D., Luhman, K, & Liu, M. 2007, ArXiv Astrophysics e-prints 
Apai, D., Pascucci, I., Bouwman, J., et al. 2005, Science, 310, 834 
Appenzeller, I., Fricke, K, Furtig, W., et al. 1998, The Messenger, 94, 1 
Baraffe, I., Chabrier, G, Allard, F, & Hauschildt, P. H. 2002, A&A, 382, 563 
Barrado y Navascues, D., Mohanty, S., & Jayawardhana, R. 2004, ApJ, 604, 284 
Basri, G. & Reiners, A. 2006, AJ, 132, 663 

Bockelee-Morvan, D., Gautier, D., Hersant. F, Hure, J.-M., & Robert, F. 2002, A&A, 384, 1107 

Bouwman, J., Meeus, G, de Koter, A., et al. 2001, A&A, 375, 950 

Brandner, W., Rousset, G, Lenzen, R., et al. 2002, The Messenger, 107 

Cabrit, S., Edwards, S., Strom, S. E., & Strom, K. M. 1990, ApJ, 354, 687 

Calvet, N., D'Alessio, P., Watson, D. M., et al. 2005, ApJ, 630, L185 

Chiang, E. I., Joung, M. K, Creech-Eakman, M. J., et al. 2001, ApJ, 547, 1077 

Ciesla, F. J. 2007, Science, 318, 613 

Comeron, F, Fernandez, M., Baraffe, I., Neuhauser, R., & Kaas, A. A. 2003, A&A, 406, 1001 
Devillard, N. 1997, The Messenger, 87, 19 

Dorschner, J., Begemann, B., Henning, T, Jaeger, C, & Mutschke, H. 1995, A&A, 300, 503 
Dullemond, C. P., Apai, D., & Walch, S. 2006, ApJ, 640, L67 
Fernandez, M. & Comeron, F. 2001, A&A, 380, 264 
Fernandez, M. & Comeron, F. 2005, A&A, 440, 1119 



H. Bouy et al.: Structural and compositional properties of brown dwarf disks: the case of 2MASS J04442713+2512164 



Table 2. Photometry 





Flux 


Ref. 
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1.662 
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Table 3. 


Pseudo-equivalent widths of emission lines (in A) 


Line 


Keck/HIRES 


CAHA/TWIN 




(R=31000, 2006 Oct. 13) 


(R=3 300, 2006 Nov. 11) 


O i 6300 


48.3±4.5 


72+3 


O i 6364 


15.8+ 0.5 


17±2 


Ha 6563 


620±330> 


124±20 


Nn6581 


5.2+1.4 


8.7+2.0 


Sn6717 


13.0±0.6 


14±3 


Sn6731 


22.5±1.4 


30±2 


Can 8498 


5.4+0.3 




Can 8542 


3.4±0.1 




Can 8662 


2.6±0.2 





'The equivalent width of such a strong, broad and asymetric line is 
only indicative 



Table 4. VLT/VISIR Filters properties 



Filter 


Central Wavelength 
Lum] 


Half-band width 
[}im\ 


PAH1 


8.59 


0.42 


Arm 


8.99 


0.14 


Siv 


10.49 


0.16 


PAH2 


11.25 


0.59 


Nen 


12.81 


0.21 


Qi 


17.65 


0.83 


Q2 


18.72 


0.88 
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Fig. 1 Emission lines in the Keck/HIRES spectra. Most of the lines are asymetric. 
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X [iim] 

Fig. 2 Low resolution CAHA/TWIN optical spectrum of 2MASS J044427 13+25 12164 (black), compared to vB8 (M7, light blue) 
and vBlO (M8, red) field dwarfs observed with the ALFOSC spectrograph on the Northern Optical Telescope at a similar resolution. 
2MASS J044427 13+25 12164 must have a spectral type intermediate between these two references. 
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Fig. 3 VLT/FORS2 image obtained in the S n filter. The target, indicated by a square, is unresolved in this 2 min exposure. The zoom 
around the target with contour plots shows that there is no clear evidence of a jet. The average seeing in this image is 1". The scale 
and orientation are indicated. 
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Fig. 4 Spectral Energy Distribution of the target (black) and best fit MCFOST model including the treatment of silicates (blue curve) 
and without silicates (red curve). The inset shows a zoom on the IRS spectrum. Some crystalline silicate features are indicated. 
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Table 5. Properties of 2MASS J044427 13+25 12164 



Stellar properties 



SpT 


M7.25±0.25 


v sin i 


12±2kms~' 


veiling r 6V5oA 


0.4±0.1 


A v 


0.0 mag 


logM [M Q /yr] 


-9.7 



Table 6. Mineralogical composition of the inner disk probed by our analysis 
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Table 7. Fractional masses of the two components 
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Fig. 5 Synthetic SED computed with MCFOST for R out =300, 10, 1 AU (blue, red, and green, respectively), all the other parameters 
being free. Fits of similar goodness can be obtained for any values of R out >10 AU. When the outer radius becomes too small, the 
disk becomes optically thick and the flux in the sub-mm/mm range decreases dramatically. 
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Fig. 6 Distribution of probabilities of the output parameters of the MCFOST fit. The inclination,/?, and R; n have the best constraints, 
the other parameters being only loosely or not constrained. 
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Fig. 7 Upper Panels: Spitzer IRS spectrum (black), and the assumed continuum estimated using the best MCFOST fit (left panel) 
and a polynomial fit (right panel) over-plotted in red, and the continuum subtracted flux (blue). Lower panels: continuum subtracted 
flux (black) and compositional analysis of the silicate emission. The blue line corresponds to the high temperature component, the 
red line to the low temperature component, and the green line to the sum of the two. 
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Fig. 8 Cumulative flux at 10 fim (black), 20 fj.m (red), 35 fim (blue) 350 //m (green), 450 im \ (cyan) and 1300 fim (magenta) as a 
function of the radius. Fluxes are calculated using the MCFOST models of iPinte et all (120061) for 2MASS J044427 13+25 12 164 best 
fit model. Of the different fluxes, 90% (horizontal dashed line) comes from within AU at 10 fim, 1~2 AU at 10 and 20 //m, and 
>20 AU for >350 fxm. 
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Fig. 9 Grain evolution on stellar p arameter: Cryst alline mass fraction of the inner disk (<2 AU) as a function of the mass for stellar 
and substellar objects (Fig. from Apai et al. 2005, updated with new measurements). 2MASS J04442713+2512164is represented 
with a red flower. 
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Fig. 10 Grain evolution on disk parameter: crystal line mass fraction of the inner disk (<2 AU) as a function of t he accretion luminos- 
ity (L acc oc x M) for Herbig Ae/Be stars (stars Ivan Boekel et al.l l2005) and brown dwarfs dApai et ail|2005l filled circle if a mass 
accretion rate has been measured, and open right-triangle if only a limit is available) and TTauri stars ( Sarge nt et alj|2006l) for the 
objects represented in Fig. [9] 2MASS J044427 13+25 12 164 is represented with a red flower. In the current state of the observations 
there seem to be no correlation b etween the mass accretion r ate and the crystallinit y on either side of the substell a r limit. Mass 
accret i on rate measuremen ts from iGarcia Lopez et all d2006l) : iMohantv etal] d2005l) : iMuzerolle et all (|2005, 2000); iGuedel et all 
(l2007h : ICalvet et all (l2005h and this work. 



